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We discovered numerically a scaling law obeyed by the amplitude of collective motion in large 
populations of chaotic elements. Our analysis strongly suggests that such populations generi- 
cally exhibit collective motion in the presence of interaction, however weak it may be. A phase 
diagram for the collective motion, which is characterized by peculiar structures similar to Arnold 
tongues, is obtained. 



I. INTRODUCTION 

Systems consisting of large populations of interact- 
ing dynamical elements are widely distributed in nature 
from communities of ants to biological cell assemblies and 
neural networks. Many such populations are similar in 
the respect that they exhibit various kinds of collective 
behavior. One possible approach to the mathematical 
study of such collective behavior is to concentrate on cer- 
tain idealized models such as interacting limit cycles and 
chaotic maps. An abundance of literature has been de- 
voted to the study of these models [@PPflfl||j7|J|,^l ■ It is 
suggested in some foregoing studies jMjTpMJ] that a pop- 
ulation of this type as a whole exhibits low dimensional 
behavior. This seems to be true even when the individual 
elements appears to be mutually uncorrelated, a situation 
in which it would seem that the population could only be 
regarded as a dynamical system of extremely high dimen- 
sion. It still remains unclear whether such collective be- 
havior could be understood in terms of low-dimensional 
dynamical systems. 

Without directly attacking the above stated problem, 
we take in the present Paper a slightly different approach 
in which we focus on a non-specific property of collective 
motion, a long-time average of the amplitude of collective 
motion. As we find, this approach, with which we in- 
tentionally avoid addressing specific features of collective 
motion, turns out to be extremely useful. As model sys- 
tems, we have chosen globally coupled tent maps, because 
these systems are particularly well suited for detailed nu- 
merical analysis. Using these maps, we have discovered 
a scaling law characterizing the amplitude of collective 
motion which holds for a particular series of parameter 
values. Further, we argue on the basis of numerical ev- 
idence that this type of collective motion should occur 
generically if an interaction exists, no matter how weak 
this interaction may be. Finally, a phase diagram for 
the collective motion is obtained, and the existence in 



it of peculiar structures similar to Arnold tongues with 
various scales is confirmed. 



II. MODEL : GLOBALLY COUPLED TENT MAPS 

Globally coupled maps (GCM) are given by an assem- 
bly of N elements whose behavior is determined by N 
identical maps with all-to-all coupling. The individual 
elements are then under the influence of a common in- 
ternal field which may be referred to as a mean field 
We assume that a single isolated element evolves accord- 
ing to X n+ i = f(X n ), where n designates discrete time 
steps. Under the interaction through the mean field h n , 
the i-th element is then assumed to evolve as: 



(1) 



where K is the coupling strength. In this Paper, we 
consider the situation in which / is a tent map: 



f(X) = -a\X\ 



a - 1 



and the mean field h n is defined as: 



1 N 



(2) 



(3) 



Thus, our system is characterized by the two parameters 
a and K in addition to the total number of elements 
N . Each tent map has a band splitting point a = \J~2 
and we thus stipulate that a satisfies \/2 < a < 2. It 
is thereby ensured that the population will never split 
into sub-populations. If the system size A is a finite, 
this finiteness becomes the source of fluctuations of the 
mean field. Such an effect may obscure pure collective 
motion. Thus we work with the limit of large N. In fact, 
we confirmed that finite size effects can be regarded as 
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noise acting on the pure collective motion. For ./V — > oo, 
the population dynamics of GCM can be described by 
the Frobenius-Perron equation ]Tl| ] for the distribution 
p{X): 

Pn+i(X) = Js(X - f[X') - Kh n )p n (X')dX', (4) 

h n = J f(X')p n (X')dX'. (5) 

We worked out a numerical scheme for the exact inte- 
gration of Eqs.(|J) and (|j) whose precision is limited only 
by round-off errors. Our scheme is almost the same as 
that proposed by Morita (To), although these numerical 
schemes were developed independently. We have con- 
firmed that collective motion appearing in the system 
described by Eqs.(Q) and (||) is almost independent of 
the initial distribution. For this reason we present nu- 
merical results in this Paper were obtained using only 
a single initial distribution, uniform over the interval 
[/°/(0):/(0)]. 



III. COLLECTIVE MOTION IN GLOBALLY 

COUPLED TENT MAPS 

In some previous works it has been argued that 

the model constituted by Eqs^j) and (||) does not exhibit 
collective motion. However, our careful analysis yields a 
contradictory conclusion. Collective motion can be ob- 
served through the dynamics of the order parameter h n 
given by Eq.(|5|). Figures [l](a) and 0(b) are return maps 
of h n . There we see that the fluctuation of h n under- 
goes quasi-periodic motion for small K, but for larger K 
it displays more complicated motion possibly with fine 
structure. Roughly speaking, larger K values results in 
more complicated collective behavior. 

The observed collective motion is similar to that found 
in other models |||7|]|] , except that the scale of the collec- 
tive motion for the present model is much smaller. Be- 
cause of this smallness, the previous works @,|l^] have 
failed to detect the collective motion for the present 
model. Based on their conclusion that collective motion 
does not exist here, they conjectured that the presence 
of window structures for the elements (which tent maps 
never possess) are necessary for the occurrence of collec- 
tive motion, but our results shows this is not the case. 

IV. PARAMETER DEPENDENCE OF THE 
COLLECTIVE MOTION : HILLY STRUCTURE 

We concentrate on the amplitude F of the collective 
motion without going into any detailed structure of the 
dynamics. We have 

F = V< {K- < h n >)2 >, (6) 



where <> represents a long-time average. In Fig.g, the 
dependence of F on a is shown for two values of K. The 
horizontal axis represents m = — log Q (2 — a), which is a 
monotonically increasing function of a on (0,2), tending 
to oo as a — ► 2. The figure shows the range of many 
hills of F, which we call 'hilly structure'. As the value of 
K is decreased, the number of the hills becomes larger, 
while the average size of each hill becomes smaller. For 
larger a (i.e. larger m), the hilly structure is difficult to 
discern in this figure, but its persistence can be confirmed 
by magnifying the scale. Hilly structure seems to exist 
for arbitrarily small values of K. Furthermore, based 
on this figure we expect that for a given value of K, 
collective motion exists for almost all a. These points 
will be discussed in further detail below. 



V. GOLDEN AND SILVER VALUES OF 
PARAMETER "a" 

We find that the role of a particular series of parameter 
values of a which we call 'golden values', play a cru- 
cial role in the formation of such hilly structure. We find 
in Section VII that these parameter values are situated in 
the middle of each hill for sufficiently small values of K. 
Some examples are indicated in Fig.^| with the notation 
.93 ~ 37- When a is identical to one of the golden values, 
one isolated tent map possesses the property that a tra- 
jectory beginning at the peak C of the tent map returns 
to C after p steps. We denote a golden value defined 
in this way as g p . Incidentally, we use the term 'golden 
value' because the value <?3 equals the golden mean. 

We refer to another special series of parameter values 
of a as 'silver values'. These are associated with the val- 
leys between neighboring hills, as argued in Section VIII. 
For instance, the silver values S4 ~ s-j are indicated in 
Fig.^. When a is identical to s p , a trajectory beginning 
at C falls into a fixed point of the tent map after p steps. 
Golden and silver values exist densely if we allow p to 
take all natural number values, while they occupy only a 
vanishing measure on the line of a. We note that g p and 
s p are generally not unique because of the multiplicity of 
p-periodic orbits, whose number increases almost expo- 
nentially with p. Despite this non- uniqueness, g p {s p ) is 
used below to represent a single golden (silver) value. 

VI. SCALING LAW OF COLLECTIVE MOTION 
AT GOLDEN VALUES 

We now consider the collective motion as a function 
of the coupling strength K, with a confined to golden 
values. We found numerically that the amplitude F of 
the collective motion obeys the law 

KF-e-Tt, (7) 
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where a is a positive constant depending on the golden 
values in question. For any golden value, the above scal- 
ing form seems to hold, as long as K is sufficiently small. 
In FigJ|, the relation between KF and K~ l is displayed 
for those golden values which are situated in the mid- 
dle of some representative hills in Fig.||. The formula in 
Eq.(^) reveals some important properties of the collec- 
tive motion. First, the collective motion persists even 
as K — > 0. Second, the observed quasi-periodic motion 
does not appear through the conventional route of Hopf 
bifurcation. This follows from the fact that if the quasi- 
periodicity were due to a Hopf bifurcation, then we would 
find that F ~ \J\K — K c \ above a bifurcation point K c . 
This is clearly in contradiction with Eq. (f7j) . In our study 
we could not even find an indication of a bifurcation. 
Third, the simple relation in Eq.(Q) survives even when 
the collective motion becomes more complicated. Note 
that different modes of collective motion as illustrated in 
Figs.|l|(a) and (b), lie on a common line of Fig.||. 

VII. FORMATION OF HILLY STRUCTURE 

The appearance of hilly structure, seen in Fig.|^, can 
be understood from two arguments, concerning hills and 
valleys. The following discussion regarding hills based on 
Eq.(Q) partially explains how collective motion comes to 
characterize almost all values of a, not just the golden 
values. 

At a given value of K each hill of F contains a number 
of golden values. One of these golden values is found to 
be representative of this hill. We call this a 'key golden 
value'. In order to define key golden values, we studied 
the growth of the width of hills with the decrease of K. 
Figure [I] shows how the width W of the hill around a 
given golden value changes with K . The width W 
around a given g p is defined by W = |at — or|, where 
a,L and an are the two values of a on either side of g p 
such that the amplitude F at a — and a = cir is 
half of F(g p ). The golden values chosen here are the 
same ones as in Figs.[| and |3[ In this figure, we see that 
W behaves as a function of K similar to that in Eq. ([?]). 
Actually, the width W for a given g p is proportional to 
KF calculated for the same g p , as long as K is sufficiently 
small. This implies that the hill narrows down to a point, 
corresponding to a = g p . Based on this relation, we 
define 'key golden values' of a given hill which exists at a 
given Kq as follows: With the decrease of K from Kq, the 
width of a hill narrows down to a single point according 
to the growth law for W. The value of a for this single 
point is that of the key golden value of this hill. The 
definition of key golden values will be refined as golden 
values displaying the scaling law Eq.(^) for values of K 
smaller than Kq. 

The formation of the hills may roughly be understood 
from the following argument. At a golden value g p , the 
invariant measure for a single map, i.e. p(X) at the equi- 



librium solution of the Frobenius-Perron equation with 
K = 0, has the rather simple form jl2) of a (p — 2)-step 
function, as displayed in FigJ|(a). As a deviates slightly 
from g pi the change of the invariant measure remains 
small, although there may appear some fine structure. 
When a weak interaction is switched on, the solution 
p n (X) of Eqs.(Q) and (||) with a — g p continues to still as- 
sume a shape quite similar to the invariant measure p(X) 
mentioned above, but it is accompanied by intermittent 
peaks (see Fig.||(b)) whose widths are roughly given by 
KF. The situation remains almost the same when a 
differs slightly from g p . If the effect of these peaks on 
the amplitude of the collective motion is more important 
than the change in the invariant measure caused by a 
small deviation of a from g Pl the statistical average F 
will be insensitive to the value of a around g p . If the 
deviation of a from g p becomes too large, however, the 
corresponding F will be dominated by the influence of 
other golden values corresponding to other hills. 

Hills formed in the manner described above are consis- 
tent with hills obeying the growth relation of the width 
W (see FigQ. Further, the above reasoning applies to 
every golden value. Thus every golden value can be con- 
sidered as a key golden value around which its own hill 
is formed. This conclusion may appear to be inconsis- 
tent with the fact that in a given hill there are an infinite 
number of golden values. Actually, for a given K, most 
golden values in a given hill are not key golden values. 
We have found that around such a K value, the ampli- 
tude F for these golden values does not obey the scaling 
law of F. This will be discussed in Section IX. It is now 
clear that collective motion is not confined only to values 
corresponding to golden values. It occurs over the inter- 
vals corresponding to hills existing around these golden 
values. 



VIII. COLLECTIVE MOTION IS DOMINANT : 
PROPERTY AROUND VALLEYS 

The hilly structure seen in Fig|| can be understood 
by analyzing the structure of F in the valleys between 
two neighboring hills. From this analysis we are able 
to conclude that collective motion occurs generically in 
the a-K plane. We specifically investigate the parameter 
dependence of F in valleys. As is seen from FigJ|, we 
find that some silver values lie precisely at the minimum 
points of F in the valleys. Figure || shows that when the 
value of a is near a silver value s p , we obtain the following 
relation for sufficiently small values of K: 

F(x\a-s p f ((3~1). (8) 

It is clear from this relation that collective motion dis- 
appears only at the point a = s p Jl3| . This relation 
seems to hold for any silver value when the value of K is 
sufficiently small. This fact strongly suggest that the to- 
tal measure of values of a at which collective motion docs 
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not exists is vanishingly small, because the total measure 
corresponding to silver values is vanishing, although they 
exist densely on the line of a. The remaining values of a 
are characterized by hills which extend from key golden 
values with widths which increase with K. Hence the 
occurrence of the collective motion is dominant over the 
parameter space of a. Because of the scaling law existing 
at golden values, this situation holds even for vanishingly 
small values of K. 



IX. MERGING OF HILLS 

Figure ^ also indicates that the valleys disappear for 
sufficiently large values of K. For each silver value, we 
find a critical value K c for the disappearance of the cor- 
responding valley. When K exceeds such a K c , a valley 
which existed at a certain silver value up to K = K c 
suddenly disappears. This occurs because two hills ly- 
ing on either side of this silver value merged into one at 
K = K c . As a result, the amplitude F at this value of 
a exhibits a sudden increase. As described above, each 
hill possesses a key golden value. After the merging of 
two hills, one of the two key golden values characteriz- 
ing the previously existing hills, say g pi , ceases to be a 
key golden value, while the other persists. As a result 
the if-dependence of F at a = g Pl changes, and typi- 
cally comes to behave as illustrated in Fig|7| We see in 
this figure that the X-dependence of F(g^) is maintained, 
while that of F(gn) changes above K c , where it comes 
to exhibit if-dependence similar to F{g±). 

The merging of two hills leads to a sudden expansion 
of the width of one of the hills, accompanied by the dis- 
appearance of the other. If a hill associated with some 
key golden value g p successively absorbs the neighboring 
hills with the increase of K, then the width W of this 
hill also grows successively. On the other hand, we have 
already seen in FigJ^ a growth law for W. These obser- 
vations suggest that the values of K c characterizing the 
disappearance of valleys around silver values should be 
correlated with our growth law for W . In the case of 
Fig.||, the golden value g± persists as a key golden value 
over a wide range of K. We calculated the values of K c 
for a series of silver values. Each of these silver values 
was chosen in such a way that in a certain range of K, 
it lies at the minimum point between the hill associated 
with #4 and a neighboring hill. In fact, there should exist 
many more silver values of this kind for a given 34 than 
those displayed in this figure. We see from this figure 
that the value of K c for a given silver value depends on 
the distance of this silver value from (74, and that this 
dependence is consistent with the growth law for W seen 
in Fig.|]. In many cases, merging occurs between hills of 
vastly different sizes. In this case, the collective motion 
characterizing the smaller hill seems to be replaced sud- 
denly by that for the larger one. In such cases, the hill 
width may appear to grow almost continuously. 



From our study presented so far, the hilly structure 
of the amplitude F of the collective motion can be un- 
derstood. This structure consists of a series of hills and 
valleys. Only at a minimum point of F in each valley, 
collective motion disappears. Thus the hilly structure 
implies that under a given coupling strength K the col- 
lective motion occurs inevitably for almost all values of 
a. Furthermore, the following properties for the collec- 
tive motion are now clear: There is a scaling law for the 
amplitude F, a growth law for the hill width W, and a 
successive merging of hills resulting in a decrease in the 
number of the hills and hence the number of key golden 
values. Clearly these features are interrelated. 

X. PHASE DIAGRAM 

We now consider a phase diagram of collective motion 
in the a-K~ x plane from the viewpoint of the amplitude 
F, For each golden value g p , we define a phase as the re- 
gion of a hill (defined by a values of F greater than some 
threshold value F t h) associated with a certain g p . Here, 
the definition is meaningful only if g p represents a key 
golden value. Since golden values exist densely, such a 
diagram will be split into an infinite number of domains 
as K — > 0, each reducing to a single point corresponding 
to a golden value. The extension of each domain becomes 
larger with the increase of K, as implied from Eq. ((?]), and 
this necessarily results in the successive merging of do- 
mains, leaving only a few for sufficiently large K. Thus 
our phase diagram resembles that for the phase locking 
between oscillators. The latter is also characterized by 
the merging and splitting of domains, in this case syn- 
chronized domains, referred to as Arnold tongues. We 
note, however, there is one property of our phase dia- 
gram which distinguishes it from that in the case of phase 
locking. In the present case, in the Fth ~ * limit, the 
domains of collective motion comprise all but a measure 
set of values of a, even for vanishingly weak coupling. 

In Fig.^|, we show a blowup of the phase diagram which 
contains only those domains of collective motion for g p 
with p < 1 1 . The numerical procedure for obtaining this 
picture essentially follows that for Fig^. We determined 
the domain for each g p with a given K first, displaying it 
as a painted region of a. We then did the same for vari- 
ous K, thereby producing FigJ^. Note that the domains 
separated for a given K correspond to distinct hills with 
distinct g p . As expected, a pair of neighboring domains 
merge at a certain value of K at which a valley disap- 
pears. Thus, the number of separated domains decreases 
with the increase of K . After several merging events, we 
are left with only one painted domain associated with g±. 
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XI. DISCUSSION : ANY UNIVERSALITY 
CLASS? 

The foregoing arguments suggest that collective mo- 
tion occurs generically for arbitrarily weak coupling. 
Such a result seems to be closely related with some global 
properties of the individual maps such as the topological 
arrangement of golden values. It is our conjecture that 
the same conclusion holds for such globally coupled uni- 
modal maps satisfying \f'(X)\ > 1 for any X. Although 
not reported in the present Paper, we also confirmed that 
similar scaling law holds for varieties of mean fields h n . 
From these studies it is seen that the analytical approach 
by Ershov and Potapov j[4| gives estimated values of F 
which are much smaller than those found in the present 
numerical results. 

For the case not satisfying the above condition for 
for instance globally coupled logistic maps, the 
present argument is expected to be qualitatively applica- 
ble under some restrictions. This is mainly because here, 
golden values are located in windows, and the stability 
of a periodic orbit there causes the so-called clustering 
phenomena [|| which are out of the present arguments. 
On the other hand, collective motion resembling to the 
present one is also observed for globally coupled logistic 
maps, where the amplitude of this collective motion is 
proportional to the coupling strength . An analytical 
study 1 15 1 also supports this idea. These suggests the 
existence of a certain scaling law of F which implies the 
qualitative applicability of our argument. 

Collective motion with a quasi-periodic property is also 
observed in spatially extended systems ||. Since GCM 
is an idealized model of a long-range coupled system, it 
does not take into account anything resembling spatial 
extension. We wish to be able to find a way to include 
the effects of spatial degrees of freedom in the descrip- 
tion of collective motion. It is hoped that the discovery 
presented here will be relevant to such spatially extended 
systems also. 

We have concentrated on the amplitude F of the col- 
lective motion without going into its phase-space struc- 
ture. The scaling relation discovered for this non-specific 
quantity holds over a wide range of parameter values for 
which the type of collective motion changes in various 
ways. Based on this point, we suspect that the collec- 
tive motion discussed in the present Paper might pos- 
sess quite unusual properties not shared by conventional 
low-dimensional dynamical systems. Detailed numerical 
study of the nature of various types of collective motion, 
which is now in progress, will clarify this point. 
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cussion, and to Y.Kuramoto, B.Hinrichs and K.Kaneko 
for encouragement and critical reading. This research is 
supported by JSPS Research Fellowships. 
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FIG. 1. Two types of collective motion at a — g$ are dis- 
played by the return map of h n . (a) For K — 0.1, we find a 
torus representing quasi-periodic motion, (b) For K — 1.0, 
the return map indicates more complicated torus-like motion 
possibly accompanied by fine structure. 



FIG. 2. Amplitude of collective motion F vs m for two 
values of K. Here m = — log a (2 — a). There are several 
hills and valleys, whose numbers increases as K decreases. 
Some golden and silver values, denoted by g p and s p , are also 
displayed for reference. 



FIG. 3. Scaling of F with K. Here KF vs l/K is shown 
for five golden values Op, p — 3 ~ 7. The values of g p are the 
same as those in FigM The data clearly exhibits the linear 
dependence of log(_K\F) on l/K. 



FIG. 4. Widths of the hills around a few key golden values. 
Here (hill width) vs l/K is shown for five golden values g p 
(p = 3 ~ 7) whose values are the same as those in Fig|| For 
each golden value, we find a good correspondence between the 
line here and that in Fig.H 



FIG. 5. (a) Invariant measure p(X) and (b) a snapshot of 
the distribution p n (X) at a certain time step n, where a is 
set at the golden value g$. (a) The invariant measure is given 
by a 3-step function, (b) The instantaneous shape of the 
distribution for K = 0.1 is similar to the invariant measure, 
except that it is accompanied by peaks with finite widths 
which appear near each corner of the steps. The snapshots 
of p n (X) for all values of a in the neighborhood of gs have a 
similar property. 
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FIG. 6. Amplitude of collective motion F near a silver 
value S7 is displayed on a logarithmic scale, where Aa = a—sr- 
The picture on the left is for values of a smaller than S7, 
while that on the right is for a larger than S7. For smaller K, 
F oc |Aa|, so that F decays to zero as a — > S7. On the other 
hand, F remains finite for larger K over the range of a in this 
figure. 



FIG. 7. KF vs K~ x for two golden values whose values of 
a differ only slightly. The plotting method is the same as that 
for Fig.^. The value of 54 is the same as that given in Fig.^. 
Each golden value is regarded as key golden value while K is 
sufficiently small. However, only the collective motion at one 
golden value, #4, is characterized by the same scaling law in 
both the small- K and large- A' regions. The K~ x dependence 
of KF at the other golden value, gn, changes and this depen- 
dence comes to resemble that at 54. This implies that the hill 
associated with 34 absorbs the hill associated with g\i. Thus, 
the golden value g4 continues to be a key golden value, while 
gn cease to be such. 

FIG. 8. The point K c at which the minimum point in the 
valley closest to the key golden value <?4 vanishes is investi- 
gated. Here Aa = a — (74. The value of (74 is the same as that 
given in Fig.^. Note that the horizontal scale is logarithmic. 
The picture on the left is for the values of a smaller than 34, 
while that on the right is for larger values. When a minimum 
situated at the silver value disappears at some K = K c , the 
width of the hill associated with 34 becomes expanded. The 
values of K c for those valleys which disappear successively 
show a logarithmic dependence on Aa. Although this figure 
ignores many other valleys which are actually present, we be- 
lieve that the manner in which valleys disappear successively 
is consistent with the variation of the hill width with K shown 
in Figg 

FIG. 9. A blowup of the phase diagram for the collective 
motion in m-K~ x space, where m = — log a (2 — a) (restricted 
in 2.82 < m < 3.22). The value of the 34 in this figure is the 
same as that in Fig.^. Each domain filled with black has been 
determined from the calculation of the hill width associated 
with the corresponding golden value. The golden values rep- 
resented here are g p with p < 11. The neighboring domains 
overlap for larger values oi K. A precise phase diagram split 
into an infinite number of domains could be imagined as an 
extension of this figure. 
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